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ABSTRACT 



Context. Hydrodynamical simulations indicate that substantial fraction of baryons in the Universe remains in a diffuse component - Warm-Hot 
Intergalactic Medium (WHIM). 

Aims. To determine physical properties (spatial distribution, temperature and density) of WHIM. 

Mefhods. Spatial structure of the soft extended X-ray emission surrounding field galaxies is carefully investigated using the XMM-Newton 
EPIC/MOS observations. Angular correlations between the galaxy distribution and the soft X-ray background extending over several arcmin 
result from the large scale clustering of galaxies. At small angular scales (below ~ 2') the excess of the soft flux is interpreted as the genuine 
emission from halos of the WHIM surrounding individual galaxies. 

Results. Bulk parameters of the average WHIM halos associated with galaxies in the sample are estimated. Typical halo has a characteristic 
radius of ~ 0.3 Mpc and a mass of 4 - 7 • 10" Mq. The average density of the WHIM in the local universe amounts to 7 - 1 1 • 10"^^ g cm"^ 
(nwHiM = 0.7-1.2%). 

Conclusions. Observations of the X-ray WHIM emission are in good agreement with the numerical simulations, but accuracy of the observa- 
tional material is insufficient to constrain the theory of WHIM. A series of deep observations of a moderately numerous sample of low redshift 
galaxies with high resolution instruments of Chandra would significantly improve our estimates of the WHIM parameters. 

Key words. X-rays: diffuse background - intergalactic medium - X-rays: galaxies 



1. Introduction 

Both theoretical arguments and observational data point out 
that a major fraction of baryons in the local universe remains in 
the form of the diffuse component. Following the first hydro- 
dynamical simulations bv Cen & Ostrike r ( 1999). several in- 
dependent groups (e.g. lDave et al. 11200 itlBrvan & VoitfeoOlt 
ICrof t et al. 2001 ) investigated the evolution of the primordial 
gas distributed outside large concentrations of mass, viz. galax- 
ies and clusters. According to simulations 30-40 % of the bary- 
onic matter has not yet concentrated in gravitationally bound 
objects. This diffuse component gradually flows toward poten- 
tial wells created by (mostly) non-baryonic dark matter. 

A rate at which the gas flows toward the galaxies depends 
on the structure of the gravitational potential and the efficiency 
of the non-gravitati onal heating, the so-called feedback (e.g. 
IZhang & Pen '2002). This comprehensive term stands for com- 
plex processes associated with transfer of energy and matter 
from galaxies back to the intergalactic medium. Although de- 
tails of these processes are not fully recognized and modeled, 
basic conclusions drawn from the calculations by all the re- 
search groups agree. It is found that the interaction between the 
infalling gas and the galaxy generates shocks which heat the in- 
tergalactic material. In effect galaxies are suiTounded by halos 
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of plasma with temperatures and densities distinctly higher that 
in the areas not affected by the feedback. 

Simulations predict a wide range of the intergalactic plasma 
temperatures. In the local universe characteristic temperatures 
of the model gas accumulating around galaxies increase to 
10^ - 10^ K, while in the more remote regions gas remains at 
the temperatures of ~ lO'^K. Also the densities of the heated 
gas fluctuate and span a wide range of magnitudes. A substan- 
tial fraction of baryons is predicted to undergo rather moder- 
ate compressi on with a rnediaii overdensity of ~ 20 over the 
cosmic mean (iDave et al. l200lh . However, most of the X-ray- 
emitting gas has densities 100 - 1000 times greater than the 
mean ( Croft etal. 2001) . In this way, the cosmological simu- 
lations of the matter evolution indicate that a measurable frac- 
tion of baryons constitutes a specific phase of matter, so-called 
Warm-hot intergalactic medium - WHIM ( Dave et al. 20011) . 

A cloud of the WHIM might reveal its existence either 
through the intrinsic emission of the thermal radiation or 
through absorption lines induced in spectra of the background 
objects. Potentially both effects can be used to evaluate phys- 
ical parameters of the gas (spatial structure, temperature, den- 
sity) and eventually to constrain models of the WHIM evolu- 
tion. Estimates based on the simulations indicate that expected 
equivalent width of absorption lines as well as luminosities 
generated by the individual WHIM clouds are extremely low. 
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Nevertheless, detections of absorption lines produced by the 
inter galactic hot plasma have been reported by several groups 
(e.g. iTripp et all l200(i iNicastro et al. I l2005l and references 
therein). Observations of O^"^ Kor absorption lines in the spec- 
trum of the blazar MKN 421 bv iNicastro et al. I (|2005) have 
been used to make a statement on statistical characteristics of 
the WHIM and its contribution to the local cosmological mass 
density. 

Extragala c tic th ermal radiation has been searched by 
iKuntz et al. I ( 1200 11) . These authors investigated the surface 
brightness of the soft XRB in the ROSAT All-Sky Survey 
(RASS) and detected a smooth thermal component with kT - 
0.23 keV which most likely was a mixture of WHIM and the 
Galactic halo emission. Existence of gaseous halos of a radius 
30 - 40 kpc surroun ding the spiral galaxv N GC 5746 has been 
reported recently by P edersen et al. I iI2005b . The authors argue 
that the observed halo is formed from the in-flowing gas rather 
than the matter expelled in the supernova explosions. 

In the present paper we continue our WHIM studies by in- 
vestigating the soft X-ray emission generated by this phase of 
the baryonic matter The main objective of the paper is to de- 
termine physical characteristics of the emitting plasma and - 
eventually - to assess contribution of the WHIM to the to- 
tal mean space density of the baryonic matter In a series of 
earlier papers we have analyzed the structure of the soft X- 
ray background (XRB). Using the RASS maps we have de- 
tected so ft enhancements of the XRB flu x around clusters o f 
galaxies JSoltan et al. 1 1 996) and galaxies (Soltan et al. 199 r). 
Assuming a thermal origin of this excess flux, S oltan et al. I 
1I2OO2I) have shown that temperature of the emission kT ^ 
0.5 keV in agreement with the expected temperature of the 
WHIM. We have extended our analysis toward smaller angu- 
lar (and spacial) scales using the ROSAT and XMM-Newton 
pointing observations and confirmed that the surface brightness 
of the soft X-ray excess flux increases as the distance to the 
galaxv diminishes, again as expected for the WHIM structure 
llSoltanetal.l2"00i. 

In this paper spatial properties of the WHIM is investigated 
in greater details. Scrupulous analysis of the XRB structure in 
the XMM-Newton data allowed us to determine the character- 
istic size and luminosity of the X-rays emitting WHIM halo 
around field galaxies. 

Two fundamental factors complicate observations of the 
WHIM emission. First, the weak and diffuse WHIM signal is 
observed against the highly variable background produced by 
a whole variety of discrete sources. The integral background 
is dominated by point-like sources identified m ostly with 
all ki nds of AGN (e.g. Alexander et al. 2003. Worslev et al. I 
I2OO5I and references therein), leaving not much space for the 
diffuse component. To increase signal-to-noise ratio of the 
WHIM component all the detected point sources are removed 
from the data. Still, the sources below the detection threshold 
generate fluctuations which impede the WHIM search. Second, 
due to relatively low temperatures, the thermal emission by the 
WHIM is expected to be soft, typically below 1 keV, with most 
of the flux at ~ 0.5 keV. At such low energies the absorption of 
extragalactic photons by the cold gas in the Galaxy is signifi- 



cant and the search for the WHIM emission has to be limited 
to high galactic latitudes. 

In the next section a short description of the observational 
material used in the analysis is given. The method of calcula- 
tions and the raw results are presented in Sect.|3l The interpre- 
tation of the detected signal is performed in Sect.|3and|5] Brief 
discussion is given in Sect.|6l 

2. Observational material 



All the simulatioiis of th e WHIM emission (e.g. lBrvan & VoitI 
bOOlilCroft et al. I2OOII) and our earlier estimates show that the 
surface brightness of an individual WHIM cloud does not ex- 
ceed few percent of the average XRB level. To achieve mean- 
ingful signal-to-noise ratio for the WHIM emission one re- 
quires extensive set of data. An eff'ective method to measure 
the weak enhancements of the XRB around galaxies is to cal- 
culate the cross-correlation function (CCF) of the XRB and the 
sample of galaxies. To draw quantitative conclusions from the 
CCF, one has to remove from the XRB maps systematic ef- 
fects which could mimic the difiiise signal. Also the sample of 
galaxies should have well-defined statistical properties. 

2.1. EPIC MOS data 

In the present analysis the X-ray data are extracted from the 
public archive^ of the XMM-Newton EPIC/MOS observations. 
The pointin gs are selected in the same way as in the work by 
ISoltan et al. (2005) and for the full description of the data re- 
duction the reader is referred to the Sect. 2.2 of that paper. Here 
only the main points are recalled and some modifications of the 
original procedures described. 

All the data obtained wifli the MOSl and MOS2 detectors 
in a "Full Frame" mode with the thin filter have been inspected 
and only pointings "suitable" for the present investigation have 
been used. The pointings with a strong source, known extended 
source or the high particle background have been classifies as 
"unsuitable". To optimize our search of the WHIM signal, new 
(in comparison to lSoltan et al. 1(120051) ') energy bands have been 
selected: 0.3-0.5 keV,0.7-1.0keV, 1.0-1.35keV, 1.9-3.0keV 
and 3.0 -4.5 keV. Energy gaps 0.5 -0.7keV and 1.35-1. 9keV 
diminish Galactic contrib ution and the strongest in ternal fluo- 
rescent lines, respectively (Nevalainen et al. 2005). 

It was pointed bv Pradas & Kerp (2005) . that to study faint 
extended sources such as the WHIM using the XMM-Newton 
EPIC detectors "all systematic effects must be well understood 
and reliable methods to eliminate their contribution have to 
be developed". In the present approach the WHIM emission 
is analyzed taking advantage of a large number of available ob- 
servations. Several hundred pointings had been inspected and 
above 150 were accepted for further processing. A median ex- 
posure time of the observation used in the investigation is be- 
low lOks. It wa s not possible to invest igate instrumental ef- 
fects reported bv IPradas & Kerp 1 1120051) for the each pointing 
separately. However, to account for the possible irregularities 
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of the CCD performance, we h ave intro duced a correction to 
the procedure described in Sotta n et al. I (12005). Separately for 
the MOS 1 and MOS 2 cameras and for the each of our energy 
bands "a sensitivity map" has been generated. The sensitivity 
map has been constructed as follows. First, all the point-like 
sources found in the individual observation have been removed. 
Then, the counts in the detector coordinates for all the pointings 
have been added. In the same way "integral" exposure map has 
been created. Finally, the summed count distribution has been 
divided by the exposure map. Since the resultant distribution 
of countrates in pixels is based on a large number of point- 
ings, fluctuations of the cosmic signal are smoothed out. Thus, 
spatial variations of the countrates represent the mean inhomo- 
geneities of the detector sensitivity. Although this procedure is 
unable to deal with the time dependent instrumental effects, it 
allows for the coherent treatment of a large set of the EPIC 
MOS observations. 

2.2. Galaxy data 

The gal axy sample has been ext racted from the APM Galaxy 
Survey JMaddoxetal.t 199 0a. Madd oxetal.. 1990b ) using the 
NASA/IPAC Extragalactic Database (NED). All galaxies with 
magnitudes between 17 and 20 have been included into the 
investigation. This magnitude selection criterion is the same 
as in the subsample o f the APM Galaxy Survey discussed by 
iMaddox et al. 1 ( 1199 6). The present galaxy sample, drawn from 
the statistically complete survey and with well-defined mag- 
nitude limits , is m ore homogeneous than the sample used by 
ISoltan et al. I f2005) and some statistical characteristics of the 
APM galaxies determined bv .Maddox et al. .11996.) are directly 
applicable to the present investigation (see below). 

In the areas not covered by the APM Survey, galaxies ex- 
tracted from the Minnesota Automated Plate Scanner (MAPS)^ 
Catalog of the POSS I have been used. Selection criteria ap- 
plied to the the MAPS Catalogue have been constructed to gen- 
erate sample of objects with statistical characteristics resem- 
bling the APM galaxy sample. Using more than 11 00 galaxies 
common to the APM and MAPS samples, a linear transforma- 
tion between the O and E magnitudes from MAPS and bj have 
been found in the form: bj = nip -Q37(mo- itie)- Fig.^shows 
the distribution of magnitudes in both samples. Although the 
correlation indicated by the dotted line is not perfect, we expect 
that galaxies selected from the MAPS have clustering proper- 
ties and the redshift distribution similar to the APM galaxies. 

3. Correlation analysis 

The average intensity of the XRB at a distance 6 from a ran- 
domly chosen galaxy is given by the formula: 



p(0) 



2j '^pxl ■ *exp 



(1) 



where the sums extend over all pointings and all galaxies, Wcnt 
denotes the total number of counts recorded in «pxi pixels sepa- 

^ The MAPS database is supported by the University of Minnesota, 
available at_http://aps. umn.edu/ _ 
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Fig. 1. Comparison of the original bj magnitudes from the 
APM Galaxy Survey and the magnitudes derived from mo and 
niE in the MAPS Catalogue. 

rated by the angle 9 from t he galaxy, and f^xn de notes the appro- 
priate exposure time (see ISoltan et al. I (12005 ' for the detailed 
description). All the data were binned into 4" x 4" pixels. The 
total amplitude of the countrate p(ff) obtained from the set of 
actual observations is a mixture of the cosmic and local X-ray 
photons as well as charged particles. One should note, how- 
ever, that any potential systematic variations of the p with the 
separation angle 6 could result exclusively from the genuine 
changes of the average extragalactic XRB flux associated with 
the galaxy sample. Thus, the data on the absolute level of the 
contaminating counts are not crucial for the determination of 
the p(6) slope. 

The flux distribution around galaxies is related to the CCF, 
w(ff), in a standard way: 



w(e)^p{d)/p-i. 



(2) 



where p denotes the total average countrate. The average counts 
p and consequently the CCF amplitude are aff'ected by the non- 
cosmic signal. 

To isolate the XRB fraction correlated with the galaxy dis- 
tribution from the total counts and to remove residual instru- 
mental effects in the p(ff) distribution, a large number of ran- 
dom sets of "galaxy" samples have been generated. The aver- 
age Pi-and(S) countrate distribution of the randomized data pro- 
vides an estimate of the XRB component uncorrected with the 
galaxy sample. The diff'erence between both distributions: 



6p{6) = p(e) - pr,nd(0), 



(3) 



which represents the net correlated signal is shown for the soft 
(0.3 - 0.5 keV) and hard (3.0-4.5 keV) energy bands in Figs.|2l 
and|3] respectively. Note change of scale at the y-axis above the 
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Fig. 2. The countrate distribution in the soft energy band (0.3 - 
0.5 keV) vs. angular distance from the galaxy averaged over the 
galaxy sample. 
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Fig. 3. Same as Fig.|2lin the hard energy band (3.0 - 4.5 keV). 



dotted line in both figures. The distribution of countrates in the 
hard energy band is noticeably flat at separations above ~ 1', 
while in the soft band systematic decrease of the flux with the 
increasing distance is present over a wide range o f separations. 
Althou gh the same eff'ect is visible in Fig. 4 of ISoltan etal. I 
112005 ') ■ the Fig. IJlis shown here as it is based on the statisti- 
cally well-defined, complete sample of galaxies. The scatter of 



points in Fig.|5]is now reduced due to larger number of galax- 
ies involved in the analysis. Since the 6p{9) function is closely 
related to the average X-ray emission associated with a single 
galaxy (see below), the statistical completeness of the galaxy 
sample is crucial for the estimates of the WHIM emission. 

The average surface brightness excess of the XRB around 
the galaxy is produced by the source(s) "associated" with the 
galaxy itself and by the nonuniform distribution of neighbour- 
ing galaxies. So, the integral X-ray emission produced by stel- 
lar sources within the galaxy and generated in the extended 
WHIM halo surrounding the galaxy, is defined here as the 
source "associated" with the galaxy. Let s, denotes the total 
average flux produced by the galaxy in the i-th energy band, 
Wgg{6) - the galaxy ACF, and Ug - the average concentration of 
the sample galaxies on the sky. The excess XRB flux above the 
background level is approximately given by the formula: 



SpiiO) ^ Si ■ fiie) + Si ■ Wgg{9) ■ 



(4) 



where the fi{0) function describes the distribution of the galaxy 
emission in /-th band convolved with the instrument Point 
Spread Function (PSF). If the WHIM contribution to the SpiiO) 
function is negligible and the angular size of the galaxy X-ray 
emission is much smaller then the width of the PSF, we have: 



6pi{0) = Si ■ P,(0) + Si ■ Wggiff) ■ rig 



(5) 



where P,(0) is the PSF in the /-th band. In this case both com- 
ponents of the excess flux at the right hand side of Eq. |5] are 
easily identified in in Fig. |2l A strong peak at small angular 
separations (below ~ 1') is generated by the point-like galaxy 
emission, while at larger separations a smooth decrease of the 
Piiff) results from the declining density of neighbouring galax- 
ies. 

The X-ray emission generated within the galaxy either due 
to the nuclear activity or in the galactic X-ray sources (X-ray 
binaries, supernova remnants) is strongly concentrated and in 
the present sample of galaxies angular extent of the galaxy 
emission is smaller than width of the PSF. On the other hand, 
the genuine emission produced by the WHIM surrounding the 
galaxy at zero lag is expected to be substantially more extended 
than the width of the PSF. One should note that this component 
would be difficult to distinguish from the smoothed flux origi- 
nating in the neighbouring galaxies. 

3.1. The effective Point Spread Function 

Since the correlation flux 6pi{0) is generated by sources scat- 
tered over the entire field of view, the corresponding point 
source signal should be averaged in the same way. It is as- 
sumed that a sample of serendipitous quasars with redshifts 
greater than 0.3 in the investigated pointings generate coun- 
trates which mimic the shape of the point sources contribution 
to the 6pi(6). The QSOs have been extracted from the NED 
and the average countrate distribution around the sample ob- 
jects has been determined in the same way as for the galaxy 
sample. The correlation signal for the quasar sample, pQso{0), 
is shown in Fig.|4] The data points representing the countrates 
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Separation [arc sec] 

Fig. 4. The countrate distribution in the merged 5 energy bands 
vs. distance from a quasar averaged over the distant quasar 
sample (points with error bars). The curves represent fits of the 
King profile to the data: solid curve - 5 energy bands merged, 
dotted curve - 3 soft energy bands merged, dashed curve - 2 
hard energy bands merged. 

summed over all 5 energy bands are adequately approximated 
by the King function: 



Pqso(0) = Po 



1+ - 

00 



+ Pb, 



(6) 



with four simultaneously fitted parameters po, 0o, y and pi,. 
Although two shape parameters of the PSF, viz. Oq and y, de- 
pend on the energy, the effect is neglected in the present anal- 
ysis. This is because the PSF width variations over our energy 
bands are minute and do not affect noticeably the present cal- 
culations. It is demonstrated in Fig.|4] where three sets of the 
shape parameters are used to fit the observed countrate distri- 
bution. The solid curve represents the best fit of the pqso(6') to 
the sum of all 5 bands. The shape parameters of the dotted and 
dashed curves have been determined separately using coun- 
trates in three soft and two hard energy bands, respectively. The 
profiles shown in Fig. 0]have been obtained by fitting two re- 
maining parameters, po and pi,. Shapes of the PSF fitted to the 
soft and hard bands are practically indistinguishable and in the 
subsequent calculations we have used the PSF shape parame- 
ters represented by the solid curve. 

3.2. The surface brightness distribution 

At separations comparable to the width of the PSF, i.e. below 
~ 1', the countrates generated in the Spi{6) by the extended 
emission are superimposed on the relatively strong signal pro- 
duced by the point-like source identified with the galaxy. 



keV 
0.30 - 0.50 
0.70 - 1.00 
1.00 - 1.35 
1.90 - 3.00 
3.00 - 4.50 




1 10 

Separation [arc min] 



Fig. 5. The average excess countrate distribution around galax- 
ies in five energy bands. For separations above 2' the power 
law with a slope of -0.7 is assumed and only the normalization 
has been fitted to the data. At separations below 2'7 individ- 
ual points with the error bars or upper limits are shown. Some 
points have been slightly displaced for clarity. 

To subtract the point source contribution from the (5p,(0) 
countrates around galaxies, the PSF has been normalized to the 
observed distribution at separations below 10" (first three data 
points in Fig.]?}. The residual flux between (f.2 and 2°7 is plot- 
ted in Fig. |5l The 0.3 - 0.5 keV band is shown with crosses, 
0.7 -l.OkeV- triangles, 1.0- 1.35keV- squares, 1.9-3.0keV 
- full dots, and 3.0 - 4.5 keV - with stars. Error bars and upper 
limits represent statistical uncertainties at Icr level. 

At large angular separation - beyond the WHIM halo sur- 
rounding individual galaxy, enhanced emission results purely 
from the clustering of galaxies. According to Eq.0] amplitude 
of the correlated signal is defined by the galaxy ACF Although 
the angular extension of a typical halo around galaxy from 
the sample is a priori unknown, simulations indicate that the 
WH IM signal is negligible at separations a bove a few arcmin 
(e.g. Brvan & V oit 200r. lCroft et a l. '2001). Consequendy we 
tentatively assume that above 2' the amplitude of 6pi{6) is dom- 
inated by the galaxy correlation term. The ACF of th e APM 
galaxies has been determined bv lMaddox et al. I (119961) : 



w„J0) = Wo ■ 



Q-0.7 



(7) 



where wq - 0.50 for 6 in arcmin. Thus, within the framework 
of the present model, the count rate distributions 6pi(ff} at sep- 
arations greater than 2' has been approximated by the power 
law: 



6pi{0) = Ai ■ 



-0.7 



(8) 
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Table 1. X-ray fluxes correlated with galaxies 



i 


Energy 


Si 


[10-^cnt/s] 




[keV] 


total 


point 


extended 


1 


0.30 - 0.50 


18.3 ±2.6 


1.15 


17.2 


2 


0.70-1.00 


18.3 ±2.0 


1.41 


16.9 


3 


1.00-1.35 


9.4 ±1.7 


0.92 


8.5 


4 


1.90-3.00 


6.5 ±2.1 


0.77 


5.8 


5 


3.00-4.50 


5.1 ±2.2 


0.82 


4.3 



where the amplitude of the correlation signal is related to the 
average integral flux associated with the galaxy: 



Ai - Si Wo rig. 



(9) 



The parameters A, for 5 energy bands have been fitted to 
the observational data points for separations between 2' and 
14' and the best power law fits are shown with straight lines 
in the Fig.|5l The error bars indicate Icr uncertainties of the A, 
estimates generated only by statistical noise. Because the 6pi{9) 
functions are based on a large number of pointings, the scatter 
of points discernible in Figs.|2land|3]is adequately described 
by the Poisson statistics of the number of counts contributing 
to each data point. 

4. The WHIM structure 

The total fluxes correlated with the individual galaxy in 5 bands 
have been determined using Eq.|9land are given in Tabled The 
average density of the APM galaxies in the magnitude range of 
\1 < bj < 2Q Hg - 0.094 4- arcmin'^ was t aken f rom paper on 
the APM galaxy counts bv lMaddox et al. I (Il996l) . 

In the present model, flux per individual galaxy .?,- has been 
split into point-like and extended components. The average flux 
contained in the point source has been estimated by fitting the 
PSF to the SpiiO) at separations below 10" (i.e. first three points 
in the 6p\{6) distribution). The integrated flux of the extended 
component was obtained by subtraction of the point source 
from the integral signal. Both components of the correlated flux 
are given in columns 4 and 5 of Table [0 Although, the point- 
like source is conspicuous in the correlation distribution, 6pi(0), 
the data in Table^show that on the average the extended source 
is substantially stronger than the point-like source. 

One should note that in all the energy bands the detected 
fluxes ascribed to the individual galaxy are extremely small. 
With the mean exposure time for a single pointing of ~ 13 ks, 
the total average signal per galaxy in the soft bands exceeds 
just two counts, while the point source amounts to a small frac- 
tion of one count. Thus, the investigated eff'ects are absolutely 
below the detection threshold for a single observation and to 
achieve a reasonable S/N ratio, the investigation necessarily 
had to be based on a very extensive data. 

Since the galaxy ACF has roughly co nstant slope of -0.7 
down to at least 03 (Maddox et al. 1996), substantial flatten- 
ing of the 6pi{9) distribution below ~2' apparent in Fig.|5]indi- 
cates that a dominant fraction of the correlated flux is generated 
indeed in the extended source. Large uncertainties involved in 
the present analysis and visualized by the error bars in Fig. |5] 



limit our analysis to the semi-quantitative, gross estimates of 
the WHIM parameters. In the subsequent calculations we make 
assumptions which inevitably introduce simplifications into the 
WHIM modeling. Albeit some assumptions cannot be verified, 
the present data do not allow for a more detailed description. 

In particular, the WHIM distribution below 2'-3' is approx- 
imated by a well-defined, smooth halo centered on a galaxy. 
Consequently, the average surface brightness of the WHIM 
emission is equal to the amplitude of the 6pi{6). Such inter- 
pretation of the observations is favoured by the maps of the 
soft X-ray emissi on based on hydrodynamic simulations by 
ICroft et al. I (120011) . where the patches of the highest X-ray sur- 
face brightness form regular circular areas centered on galax- 
ies. Accordingly, the mean angular extent of the WHIM halo, 
^aio, is estimated using the relationship: 



Jo 



2ne5pi{e)d0^sf 



(10) 



where 6pi{9) is the excess count rate distribution after the re- 
moval of the central peak produced by the point-like compo- 
nent and .s'i^'" denotes the flux of the extended source. In rough 
agreement with Fig. |5] we assume that for the small separa- 
tions the excess countrate Spiiff) is flat and in the energy band 
0.3 - 0.5 keV, 5pi{0) ~ 3 • lO^'^cnts"' cm"^ Substituting the 
relevant quantities into Eq. ^| the implied mean halo radius 
6'haio - 2'8. Because of the multi-step procedure involved in the 
calculation, the formal error estimates of 6'haio are difficult to as- 
sess. It is likely that our "best estimate" of the 6pi{ff) amplitude 
at large separations has been overestimated due to some con- 
tribution of the individual halo to the power law section of the 
distribution. The lower value of i'i^'" would imply smaller size 
of the halo. In our opinion, realistic estimate of X-ray emitting 
halo size are between 2' and 2f8. Calculations of the baryonic 
mass of the halo described in the Sec. |5]below account for this 
uncertainty. 

4.1. Temperature of the WHIM 

Countrates of the extended source listed in the last column 
of Table ^ have been used to estimate the mean tempera- 
ture of the WHIM emission. We applied the following proce- 
dure. First, a grid of thermal plasma emission spectra based 
on the Raymond & Smith (1977) code has been generated 
for a wide range of temperatures and metal abundances. The 
MIDAS/EXS AS software has been used. The spectra have been 
calculated for redshift z = 0.12 (see below) and subject to the 
low energy absorption by cold gas in the Galaxy with the hy- 
drogen column density of 2.2 ■ 10^*^ cm"^ (the average value in 
the sample). 

Next, the spectra have been convolved with the effective 
area of the X-ray telescope/EPIC MOS detector system and the 
model counts in five energy bands have been compared with 
the data. The best fit with x^ - 10.1 for 3 dof was obtained 
for kT - 0.50 keV and unrealistically low metal abundances 
relative to the "cosmic abundances" ^ = Z/Z^osmic = 0.01. We 
conclude that the observed counts are inconsistent with a single 
temperature model - the distribution of countrates in 5 energy 
bands is too wide to be fitted by a simple thermal spectrum. 
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A presence of cluster galaxies in our sample provides a nat- 
ural explanation of the discrepancy between the data and the 
single temperature model. Although, our galaxy sample has 
been carefully examined from the point of view of the cluster 
contamination and all the known clusters have been removed 
both from the galaxy sample and the X-ray maps, it is evident 
that this procedure is not effective in elimination of all the clus- 
ters in the investigated area. In particular, our data unavoidably 
include some unspecified number of X-ray clusters with the 
brightest galaxies above the faint magnitude limit of nib = 20. 
Since the sensitivity of the present investigation is very high, 
even a moderate contamination of the galaxy sample with clus- 
ter members could introduce the measurable signal in all the 
energy bands. 

To subtract contribution of the serendipitous clusters to the 
extended emission, a thermal spectrum with kT - 5 ke V and 
metal abundances ^ = 0.20 has been adjusted to two hard en- 
ergy bands (1.9 - 3.0 and 3.0 - 4.5 keV). Then, the cluster 
spectrum has been subtracted from the remaining 3 soft energy 
bands and the the fitting procedure to the residual countrates 
has been repeated. Obviously, the temperatures of the new fits 
are substantially lower than those obtained for the original data. 
Fits with acceptable x^ have been obtained for temperatures 
in the range 0.20 - 0.30 keV with rather weak constraints on 
the metallicity. For temperatures close to kT = 0.20 keV high 
metallicities are favoured, while for kT ^ 0.25 low metallici- 
ties ^ ^ 0.30 give acceptable fits. 

5. Baryon densities 

Our estimates of the size and temperature of the average halo 
are now used to assess the characteristic densities and the total 
baryonic mass of the WHIM. In the calculations we assume 
that plasma fills the spherical halo with constant density and 
the temperature kT = 0.25 keV. According to Sect. 0] the halo 
radii of 2' and 2!8 are considered. 

Roughly 80 % of the APM galaxies in the mgnitude range 
n < bj < 20 have redshifts between 0.06 an d 0.22 with th e 
median of the distribution Zmed = 0.12. (Mad dox et al. Ill99q) . 
Luminosity and angular diameter distances defined by this red- 
shift are Dl = 560 Mpc and Da - 447 Mpc, respectively; angu- 
lar size Shaio - 2' corresponds to the radius twhim = 260 kpc = 
8.0 • 10^^ cm (Ho = 70km/s/Mpc, Q„, = 0.3, Qa = 0.7 is 
assumed throughout the paper). We now substitute the results 
obtained in the paper and all the relevant data into standard for- 
mulae: 



Vhalo = ^ ^ '•halo ' 
^halo - Vhalo ' Je(T) , 
^halo = 47rD^ ■ 5 halo , 



(11) 

(12) 
(13) 



where Vhaio, ^haio and 5 halo denote respectively: the halo vol- 
ume, luminosity and the corresponding flux at distance Dl 
from the halo; Je{T) is the volume emissivity of plasma at en- 
ergy E and temperature T. Using the basic formulae for the 



thermal Bremsstrahlung (e.g. iRvbicki & Lightman 1119851) we 
get for the 0.3 - 0.5 keV energy band: 



JEiT) = 6.38 ■ 10 



-25 ,2 



erg 



(14) 



where rip is the hydrogen density in cm"^. Eq.[2]is valid for 
the zero metallicity plasma with Y - 0.25, as both the simu- 
lations and our results favour low Z models. Finally, the coun- 
trate Sj"' = 17.2cnt/s is expressed in the cgs units using the 
conversion factor of 1 cnt - 0.65 -10"" erg/s calculated for the 
thermal spectrum with the kT = 0.25 keV. Combining all the 
numerical data we get the plasma density in the WHIM halo 
Phaio = 2.3 ■ 10"^** gem"-' and the baryonic halo mass Mhaio - 
6.9 ■ 10" M© for the ^haio = 2', while for 6lhaio = 2:8 we have 
respectively phaio = 3.9-10"2'*gcm"^ andMhaio = 4.2-10" M©. 



5.1. WHIM contribution to Obaryon 

Within the framework of the present model, the WHIM con- 
tribution to the Obaryon IS defined by the space density of the 
WHIM halos centered at galaxies. Because of a large radial 
density gradient, our magnitude limited sample of galaxies 
is not well-suited for such estimates. To correct for this ef- 
fect, the mean spatial density of galaxies contributing to the 
observed Spi(ff} distribution averaged over redshift should be 
weighted by the number of galaxies in the sample. Let n{z) de- 
notes the average number of galaxies in the magnitude range 
17 < i>y < 20 as a function of redshift. The effective spatial 
density (A^) of galaxies in the sample is given by the formula: 



<A^) = 



jN(z)-n(z)dz 
J n(z) dz 



(15) 



where N(z) is the number of galaxies in 1 Mpc"^ at redshift 
z within the apparent magnitude range 1 7 < bj < 20 . Using 
the luminosity function of galaxies by iFoIkes et al. I (119991) 
adapted to Ho = 70km/s/Mpc, we have determined the both 
redshift distributions expected in the APM sample, N{z) and 
niz). The latter one is shown in Fig. |5] superimposed on the 
observed histogram of 230 gal axies in the magnitud e range 
17 < bj < 20 as given by iMaddox et aTl (Il996l) . After 
substitution of these distributions into Eq. [21 we finally got 
(A^) = 2.4 • 10"' Mpc ^ and the average WHIM density 
PWHIM = Mb-Ao/{N} = 1.1 ■ 10"^'gcm"^ for 0haio = 2' and 
6.8 ■ lO^^'^gcm"-' for 6'haio - 2f8. Assuming the total baryon 
density pbaiyon = 0.04 pcr = 3.7 ■ 10"^' gcm"^ the WHIM de- 
tected in the soft X-rays contributes 18 - 31 % to the baryon 
density, or f^weiM = 0.7 - 1 .2 %. 



6. Discussion 

The present investigation is concentrated on two objectives. 
First, we demonstrate an existence of the extended soft emis- 
sion correlated with galaxies, which in a natural way is identi- 
fied with the thermal emission of the WHIM. Second, using the 
estimates of the X-ray flux correlated with galaxies, we deter- 
mine the physical parameters of the WHIM: the temperature, 
plasma density and the total mass of the emitting clouds. Very 
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Fig. 6. The model redshift distribution for galaxies in the mag- 
nitude range 17 < bj < 20 d erived using the luminosity 
function given by iFolkes et al. I lll999 f). The distribution has 
been normalized to th e histogram of 230 galaxies presented by 
fMaddoxetaLl ( ll9"96l) . 



low amplitude of the extended emission superimposed on the 
granular XRB generated by discrete sources, makes the analy- 
sis difficult and sensitive to the cosmic and instrumental effects. 

The first objective was achieved by the careful evaluation 
and extraction of several interfering signals and by using an ex- 
tensive observational material. The cross-correlation technique 
appears to be a highly effective tool in the present analysis. This 
is because it integrates the delicate X-ray glow suiTounding 
galaxies and efficiently eliminates various intruding signals. 
The large number of pointings and galaxies reduced the fluc- 
tuations virtually to the Poissonian limit defined by the number 
of galaxy-X-ray count pairs. Our estimates of the amplitude of 
the extended emission in 5 energy bands have been obtained 
under the assumption that the WHIM distribution is coiTelated 
with galaxies, without any further constraints. Strong correla- 
tion between the WHIM and the galaxy distribution is indicated 
by all the hydrodyiiamical simulations ("e.g. Dave et al. 2001 : 
iRrvan fc Voitl200lUCroft et al. l200lh . 

Detailed physical characteristics of the WHIM determined 
in the paper are not free from assumptions which - albeit plau- 
sible - could not be verified by the observations. Most of our 
estimates of the plasma parameters have been obtained un- 
der the assumption that the typical WHIM cloud has spheri- 
cal shape. The two-point correlation function gives the WHIM 
signal averaged over the azimuthal angle. Thus, the method is 
unable to identify filaments. Although most simulations predict 
that the WHIM distribution exhibits filamentary structure, it is 
likely that the spherical symmetry approximation describes the 
WHIM emission satisfactorily. This is because, the elongated 
structures "visible" in the WHIM simulations have extremely 
low density and surface brightness, while the regions with high 
surface brightnes s are roughly circular . This effect is visual- 
ized in Fig . 1 in iBrvan & Voit I 1120011) and Figs. 3 and 5 in 
ICroftetal.lll200llK 

In the calculations of the halo mass we assumed smooth 
and uniform distribution of the gas within a sphere surrounding 
the galaxy. Such approximation makes our calculations pos- 
sible but it certainly constitutes a simplification to the actual 
halo structure. Nevertheless, simulations provide some justi- 
fication for our approach. A qualitative inspection of the dif- 
fuse emission maps generated by iCroft et al. (2001i) reveals 



surprisingly simple structure at scales typical for the individual 
halo. Apparently the simulated structures seem to be fairly well 
modeled in the present paper Although "visual agreement" be- 
tween the maps and the simple halo model of the WHIM sup- 
ports our calculations, the range of the parameter uncertainties 
of the single halo is probably wider than that formally assessed 
in the paper. Consequently, our estimate of the average WHIM 
density also are subject to higher uncertainties. A long chain 
of calculations constrained by the series of assumptions make 
our error estimates disputable. On the other hand, our best es- 
timates of some parameters (i.e. the plasma temperature and 
density, contribution of the WHIM to the total baryon density) 
fit well to the coiTesponding quantities obtained in the hydro- 
dynamical simulations. 

One should note that to constrain the simulations, we need 
substantially more accurate observational material. It is doubt- 
ful that the WHIM emission at separations above several ar- 
cmin from the galaxy could be measured with substantially 
better precision using the present day instruments. However, 
at smaller separations, below 2' - 3', where the WHIM surface 
brightness amounts to several percent of the total soft XRB, 
a new deep observations of a moderately numerous sample 
of galaxies at low redshifts would significantly improve our 
data. Also, at separations of several arcsec one could take ad- 
vantage of the superb angular resolution of Chandra X-Ray 
Observatory and measure the WHIM emission in the imme- 
diate vicinity of galaxies. 
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